Astronomy & Astrophysics manuscript no. Galaxy5Archive 


February 1, 2008 


(DOI: will be inserted by hand later) 





Transient X-ray Emission from Normal Galactic Nuclei 
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Abstract. X-ray transients appeared in optically non-active galactic nuclei have been observed in recent years. 
The most popular model explaining this kind of phenomena is the conventional tidal disruption model. In this 
model, when a star moves within the tidal radius of a black hole, part of the star materials will fall into the black 
hole through an accretion disk, which gives rise to the luminous flare. We propose that the X-ray emission may not 
necessarily come from radiation of the accretion disk alone. Instead, it may be related to a jet. As the jet travels 
in the interstellar medium, a shock is produced and synchrotron radiation is expected. We compared the model 
light curve and the synchrotron radiation spectrum with the observed data, and find that our model explains the 
observed light curve and late-time spectrum well. Our model predicts that these transient active galactic nuclei 
could be sources of the future gamma-ray satellites, e.g. GLAST and the emission r! egion will be expanding with 
time. 
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1. Introduction 

The strong X-ray variations at the galactic nuclei of 
optically non -active galaxies has been discovered since 
1990s (e.g. seelBade et all 19961 iKomossa fc Greiner|[l999l 



iGrupe et alJ fl9991. Some of these variations are actually 
from Seyfert galaxies (e.g. WPVS0007 (|Grupe et al.lll995l 
12007ft ;) while some others are from opticall y normal, non- 
activ e galaxies (e.g. RX J1242.6-1119A) (|Halpern et al 



12004ft . In this paper we would like to focus on the emissions 
from normal galaxies. There are common characteristics 
of these sources. Firstly, all these sources have been bright 
sources, their X-ray luminosity could go up to about 10 44 
erg/s. Secondly, they have shown high level of variabil- 
ity in their X-ray light curve within years. At the 'high 
state', the luminosity of one source could be at least 100 
times higher than its luminosity at its 'low-state'. Thirdly, 
most of them have a super soft spectrum during the flare, 
with ef fective black body temperatures only abou t 10 - 
lOOeV (|Komossa fc Badelll999l Irlalpern et al.ll2004ft . The 
classic examples which satisfy these characteristics are RX 
J1624.9+7554, RX J1242.6-1119A, RX J1420+5334, RX 
J1331-3243and NGC 5905. Notejiowever, that NGC 5905 



is a starburst galaxy (jHalpern et al.ll2004l ). 

Many scenarios have been proposed to explain this 
phenomena. However, most of them fail to explain some 
of the observed results. A detailed discussion on these sce- 
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narios can be found in Komossa fc Bade (|l999ft . Among 
all the listed models, the tidal disruption model is the 
most commonly accepted model and it gives the most sat- 
isfactory explanation to the observations by considering 
the radiation from the disk. However, it still cannot fit 
the situation perfectly. 

In this paper, we try to propose an alternative radia- 
tion mechanism which may explain the areas that cannot 
be explained by the conventional tidal disruption model. 
We note that when a star is captured by a supermassive 
black hole at the galactic centre in the tidal disruption 
model, the star will be accreted to the black hole and 
a jet may be launched during the accretion process. We 
propose that the radiation emitted in a shock produced 
by the mildly relativistic jet can also form an important 
component in the X-ray flare. In section 2, we discuss the 
characteristics of some classic flares. In section 3, we will 
give a more detailed description on the tidal disruption 
model for completeness. In section 4, we provide a detailed 
discussion on our jet emission model and some numerical 
results. Finally, we give a brief conclusion and discussion 
in section 5. 

2. Typical sources 

This section summarizes the observational results of the 
commonly accepted candidates of the X-ray flare events 
satisfying the criterions given in the introduction. We fo- 
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cus on the spectral properties rather than the scale of lu- 
minosity variat ion, which has been discussed in depth a l- 
ready (See, e.g. lVaughan et al.ll2004l ; lHalpern et al.ll2004l ). 
Among all, NGC 5905, RX J1242.6-1119Ahave been stud- 
ied in detail while little information is known for RX 
J1331.9-3243. 



2.1. RX J1242.6-1119A 

This is an interesting sourc e and its detailed observa - 
tion al results can be found in iKomossa k. Greiner (1999) 



and iKomossa et all (|2004l ). RX J1242.6-1119A was ob- 
served two times by ROSAT, with the first time dur- 
ing RASS(ROSAT all-sky survey) during DEC 1990 - 
JAN 1991. The second time was a pointed observation 
with PSPC (position sensitive proportional counter). It 
was not detected in the RASS (this gives count rate 
< 0.015 cts/s) but it became visible in the pointed ob- 
servat ion with the maximum c ount rate given by 0.3 
cts/s (jKomossa fc Greinerlll999l ). From this detection, it 
is found that the photon index of this source is —5.1 ± 0.9 
if it is fitted by a power law modelled = 1-5) and gives a 
temperature of 7 x 10 5 K when fitted to a blackbody model 
(xled = 0-7) (|Komossa fc Greineilll999l ). This means that 
the source is very soft. 

There were two follow up observations by Chandra (at 
9 MAR 2001) and XMM (in 21-22 JUN 2001). The obser- 
vation by Chandra helped precise localization of the X-ray 
source. From the observation by XMM, it is found that 
the photon index of th e source has changed to — 2.5 ± 0.2 



(jKomossa et al.l 120041 ) . In other words, the source has 



hardened. Since both Chandra and XMM observations 
suggest that the detection is a point source, we believe 
that the detected X-ray was still related to the flare at 
the nucleus rather than the other parts of the galaxy 



( Komossa et al. 2004) 



Careful studies in the observation by different tele- 
scopes in the optical band suggests that th ere is no AGN 
activ ities in the galactic centre (see, e.g. lHalpern et al 



2004). There is also no radio det ection for this s ource by 
NRAO VLA Sky Survey (NVSS) ()Komossall2002h . 



2.2. NGC 5905 

NGC 5905 is not non-active since HII lines can be found 
fr om the opt i cal spectrum. However, as discussed in 



Halpern et al.l (|2004D . due to the huge variation in lumi- 



nosity, tidal disruption event is still the favoured model 
for the flare. Therefore we included this source here for 
discussion. 

NGC 5905 was observed six times by ROSAT, with 
the first three times observed during RASS, the next two 
times were pointed observations by P SPC and the last 
one by HRI (High Reso lution Imager) ( Bade et"al]|l996t 
Komossa fc Badelll999h . However, it could only be de- 
tected in one of the RASS observations, one of the pointed 
observations and the HRI observation with count rates 



found to be 0.6 cts/s, 0.007 cts/s and 0.0007 cts / s respec- 
tively (jBade et alJll996l IKomossa fc Badelll999l) . Clearly, 
there was a huge drop in luminosity. The spectrum has 
also been found to be hardened. From the observational 
result of RASS, the photon index is found to be — 4.0±0.4 
in the power law model, while in the pointed observatio n, 
the photon index become —2.4 ± 0.7 ( Bade et al. 19961) . 



NGC 5905 is observed again by Chandra in OCT 2002. 
However, the i mage is a diffused so urce, with all photons 
below 1.5keV (jHalpern et al. I l2004l) . It is impossible for 
us to determine if the photons positioned at the nuclei 
of the galaxy came from the flare or from the star burst 
region. The collected photons are also too few for spectral 
analysis. 



2.3. RX J 1624. 9+ 7554 

RX J1624.9+7554 was searched by ROSAT for two times, 
but only detected in the RASS in OCT 1990 with a count 
rate of 0.54 cts/s. The second observation was made in 
JAN 1992 by PSPC giving the co unt rate an upper limit 
of 0.023 cts/s (|Grupe et al.lll999h . Several spectral mod- 
els have been tried to fit the data but only models with 
a power law component (i.e. power law, power law free, 
power law plus a black body spectrum) can give a rea- 
sonable fit. T he photon index in these models range from 
—4.1 to —3.3 ( Grupe et al.lll999T ). In case the black body 
mo del is involved, the temperature is found to be 1.2 x 10 6 
K dGrupe et al.lll999l) . 

There was a follow-up observation for this source by 
Chandra. Since there were only 3 to 4 c ounts collected 
in th e whole detection period (~ 10 ks) ( Halpern et al.l 
120041 ). s pectral fitting i s not meaningful. However, as ar- 
gued bv lHalpern et al. ( 20041 ) . since all photons lies in the 
range of 0.7 - 4.8 keV, a blackbody spectrum is unlikely 
and so a power law spectrum can be assumed. 

Similar to RX J1242.6-1119A, studies in the observa- 
tion by different telescopes in the optical band suggest 
tha t there is no AGN ac tivities in the galactic centre (see, 
e.g. lHalpern et aTll2004h . There is also no radio detection 
for this source by NVSS (|Komossall2002lh 



2.4. RX J 1420. 4+5334 



The i nformation below is a summary of [G reiner et al 



(l2000h . This source was observed four times by ROSAT. 
However, it was detected only in the second observation, 
and the other three observations could only give an up- 
per limit to the count rates. Although there are follow up 
observations by Chandra, the results have not been pub- 
lished yet. By comparing the flux of the brightest observa- 
tion with the other three observations, this source showed 
a variation of more than 150. The spectral type of this 
source is unknown since all tested models (i.e. blackbody, 
disk blackbody and power-law) seems to fit the spectrum 
equally well. However, it is for sure that it is a soft source 
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since the photon index in the power-law fit is found to be 
-5.8. 

Note that optically there are two galaxies within the 
ROSAT error box, and we cannot single out which one is 
the host galaxy of the flare. The brighter galaxy, which is 
likely to be an elliptical or early spiral galaxy, is of redshift 
z = 0.147 ± 0.001. The fainter one, which seems to be a 
non-active galaxy, gives the lower limit on the redshift of 
z ~ 0.07. Since the optical observations suggest that its 
host galaxy is likely to be non-active, it is treated as a 
candidate for a tidal disruption ev ent. There is n o radio 
detection for this source by NVSS (|Komossall2002 ). 



sources and their spectrum can be fitted by the power law 
model. Moreover, for sources with follow-up observations, 
their spectrum seems to be hardened after the peaks of the 
flares. In general, the flares faded away in a few months 
to years. The flaring or decayi ng light curves o f the a bove 



sources can b e foun d in e.g. [V aughan et alj (|2004f ) and 



Halpern et al.l (|2004l ). lHalpern et al.l (2004 ) have shown 



that the decaying light curve can also be described by the 
power law model with the power law index given by -5/3. 
Among all the sources, NGC 5905 has the largest number 
of observations. Therefore, we will use the data from NGC 
5905 for illustration in Section 4. 



2.5. RX J 1331. 9-3243 

The information 



below is a summary of 
This source appeared two 
Its luminosity has 



Reiprich fe Greinerl (|200lD 
times in five ROSAT observations 
increased an order of magnitude within a week. Other 
details can be found in Table [1] It was also non-detectable 
six month prior and after the flare. The optical spectrum 
of this source shows no signs of AGN activities so it is 
treated as a candidate for the tidal disruption event. 

2.6. Summary of the sources 



Table 1. A summary o f the s o urces, quoted from|K omossa 
(l2002h. iGreiner et alJ (l2000h . iReipric h fc Greinerl (l200lh: 
Komossa fc Greinerl ( 19991 ) 



Galaxy 


z 


kTbb 




ix,bb 


Photon 






[kcV] 




erg/s] 


Index 


NGC 5905 


0.011 


0.06 


3 


x 10 42 


-4.0 


RX J1242.6-1119A 


0.050 


0.06 


9 


x 10 43 


-5.1 


RX J1624.9+7554 


0.064 


0.097 




- 10 44 


-3.3 


RX J1420.4+5334 


0.147* 


0.04 


8 


x 10 43 


-5.8 


RX J1331.9-3243 


0.051 






- 10 44 


-3.8 



Please refer to Section 2.4 for discussion on the redshift. 

Table [1] shows t he impor t ant p roperties of these 
sources. According to lKomossal (|2002l) . L Xj bb given in col- 
umn 4 are intrinsic luminosities of the flares in 0.1 — 2.4 
keV band based on the black body model. They should be 
a lower limit of the peak luminosities since we are most 
likely not observing the peak of the flares. The photon in- 
dexes listed in the column 5 are the one found in the 'high 
state' of the galaxy. 

As a summary, all the five sources listed above exhibit 
some common characteristics. For example, the sources 
have shown huge variation in their X-ray luminosity, but 
they are silent in optical and radio band. They are all point 



3. Review of Tidal Disruption Model 

The ideas of tidal disruption by a superm assive blac k 
hole has started as early as 1970s (e.g. see lReeslll988l) . 
It is proposed that the sudden rise in X-ray luminosity 
of the galactic nuclei of optically non- active galaxies is 
related to this type of events (e.g. see Bade et all Il996l : 
Komossa fc Greiner 19991 iGrupe et al. 19991 ). Wereview 
the classic description of the tidal disruption model very 
briefly here for completeness. 

Let Mbh be the mass of the galactic black hole, M* 
and r* be the mass an d radius of the disrupted star re- 
spectively. According to Reesl ( 19881) . if a star orbited to a 
position within the tidal radius r t of a black hole, where 



n 



5 x 10 12 



BH 



(io 



M* 



— cm, 



(1) 



the star will be heavily distorted. We consider the case 
with rt greater than the Schwarzschild radius of the 
black hole so that only part of the star will be swal- 
lowed by the black hole. Conventionally, electromagnetic 
radia tion emitted in three difference processes is consid- 
ered (jKomossal 120011 ). The processes include (1) interac- 
tion among the disr upted, but bounded stellar materi- 



als (Kim et al.l 19991 ) (2) the accretion process (e.g. see 
Komossal 120021 and the discussion below) (3) the inter- 
action between the unbounded stellar materials with the 



interstellar medium ( Khokhlov fc Melialll996t ). 

According to the numerical simulations done by 
Aval et all (|2000l ). for a solar type star, about 10% - 50% 
of the mass will get bounded and eventually accreted into 
the black hole while the remaining part of the star will be 
driven away to conserve the angular momentum. As the 
material accrete to the black hole, radiati on is given ou t 
due to the hot accretion disk. According to Ulmer ( 19991 ). 
the black body temperature of accreting a solar-typed star 
through a thick disk at the tidal radius is given by 



T,, 



2.5 x 10' 



M 



I 10 6 M(7 



BH 



1/12 



K. 



(2) 



and this will characterize the spectrum of the tidal dis- 
ruption event. However, more detailed description on the 
spectrum relies mainly on numerical simulation rather 
than analytical calculation. 
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Halpern et al.l (|2004f ) have pointed out that the decay- 
ing light curve follows the power law of the form 



L x <x 



t-t D 
tn-t 



D 



-5/3 



(3) 



where to is the time when the fla res started and tp is 
the time when the disruption began. lHalpern et al. (2004) 
have also estimated that for a tidal disruption of a non- 
rotating star, 

\ 1/2 

(to - t D ) « 1.32 x ( -^f- month. (4) 



( Mbh 
V10 6 M Q 

Th is gives the timescale for the fading of the flare. 



Halpern et al.l (J2004J) have given estimated tp using 



equation and the results are shown in Table [U The 
third column t b StP — tp gives the number of days between 
the first observation of the flare and tp. 



Secondly, accord ing to the fitting results listed in 
Halpern et al. (2004). the peak flux of two flares only oc- 



curred about tens of days before their first detections by 
ROSAT. For a flare which last over months to years, it 
would be too lucky if we could detect these sources at 
the starting time. Although one may argue that the flare 
may have a long peak duration (about 1 year) if the black 
hole is small enough so that it is easier to catch the peak 
of the flare ( Halpern et al. 20041 ). it will then be strange 
since the luminosities of all the detected sources falls in 
the follow-up detections and none stay at the peak. 

In order to explain the above points, a new radiation 
mechanism is needed. Moreover, up to now, there is no 
specific model or simulation presented to explain the ob- 
served hardening of the spectrum as time increases. We 
explore below whether the jet contribution can solve these 
problems. 



Table 2. The parameter tp of different sources 



Galaxy 


t D [yr] 


tobs.p — t D [day] 


NGC 5905 


1990.4 


« 50 


RX J1242.6-1119A 


1991.36 


« 584 


RX J1624.9+7554 


1990.7 


» 24 



In general, the luminosity-time relationship in the ob- 
se rvations can be well fitted by equation ([3]) (refer to Fig 
4, Halpern et al.ll2004[ ). The estimation of t he black hole 
mass using equation (0) is also reasonable. iDonlev et al 



(2002) lias also estimated that the frequency of the tidal 
disruption event is about 10 -5 per year in every normal 
gala xy. This frequen cy also matches with current observa- 
tion (|Komossall200ll) . 

4. Jet Emission 

4.1. Motivation 

Although the tidal disruption event explains the ob- 
served flares reasonably well, especially that it can provide 
enough energy for the huge luminosity variations and cor- 
rect decaying timescale, this pure disk emission model still 
misses several important points in observations. 

Firstly, a thermal component can usually be found in 
disk emission. However, according to past analysis of the 
sources (cf. Section 2) the thermal component can be very 
insignificant in observed spectrum. Especially, the black 
body model can only give a poor fit to RX J1624. 9+7554 
even during the flare. For the sources with follow-up detec- 
tions by either ROSAT, Chandra or XMM, it can be seen 
that generally the spectrum hardened and can only be de- 
scribed by the power law model. Therefore, it is likely that 
other radiation mechanisms may involve in the process. 



4.2. The Model 

Consider that during an accretion event, a jet can be 
formed and is ejected from the black hole/ accretion disk. 
This idea comes from the fact that accreting black holes 
are seen to be accompanied with jet. One example of such 
system is the microquasars. Studies on microquasars re- 
veal that these objects behave very differently in their 
high/low state. In their low state, there is evidence of 
jet emission although the bulk Loren tz factor of the je t 



is likely to be less than 2 (e.g. see iGallo et al.l 120031 ) 



When they are in the high/soft state, ther e is evidence 



that t he jet formation is greatly suppressed ([Fender et al. 
19991 : IGallo et al.ll2003f ). However, in their "very high" 
state, the jet reappears again. Unlike the jet seen in 
the low state, the jet is very powerful and highly rela- 
tivi stic in the " very high" state of a microquasar (e.g. 



see 



Fenderl 120031 ). This example shows that transient ac- 
creting systems may also b e accompanied with power- 
ful jets. IWang & Hul (2005) have also discussed the jet 



formation scenario which occurred during stellar capture 
events by galactic black holes to explain the X-ray cav- 
ities in some elliptical galaxies. If jet does exist dur- 
ing a tidal disruption event, the ejected jet will interact 
with the interstellar medium (ISM) and decelerate ac- 
cordingly. Synchrotron radiation from jet can explain the 
power-law decaying light curve and non-thermal spectrum 



Piranl 19991 van Paradiis et all 2000: 


Granot & Saril2002: 


Zhang k Meszarodl2004l: Zhanel 2007 


). Here we try to use 



the jet emission to explain the non-thermal X-ray flares 
from normal galactic nuclei. 

In the numerical simulation given below, we assume 
that the jet expands l aterally and satisfies four dynamical 
differential equations ( Huang et al. 2000 hi ). 



dR 
~dt 

dm 
~dk~ 



(3cr(r + Vr 2 - 1), 

; 2nR 2 (l — cos9)nm p , 



(5) 
(6) 
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dm 
dd 

~dt ' 



r 2 - 1 



M e 



2(1 - e)rm' 



R 



(7) 



(8) 



where T is the bulk Lorentz factor of the jet, R is the 
radius of the shock, j3 = v/c, where v is the bulk velocity 
and c is the speed of light in vacuum, m is the swept up 
mass by the jet, 9 is the half opening angle of the jet, n is 
the environment density, M e j is the ejected mass, e is the 
radiative efficiecy and c s is the comoving sound speed. 

Since the jet is probably mildly relativistic, as it moves 
into the ISM, a shock front is forme d and synchrotron ra- 
diatio n will be resulted. According to lRvbicki fe Lightmanl 
( 19T9h , the synchrotron radiation power at frequency v' in 
the jet comoving frame is given by 



P'iy 1 



V3e 3 B' 



7«,B 



where e is the charge of the electron, B' is the comoving 
magnetic field, m e is the electron mass, j e is the electron 
Lorentz factor. The critical frequency v' c , the minimum 
and maximum Lorentz factor 7c.mim 7e.max and the func - 
tion F(x) is given by the equations ( Huang et al.ll2000~bT ) 



7e,n 



7e,! 



Fix) 



£e (r-l)^^ + lfor P >2, 



10 8 (B'/1G) 



-1/2 



-Hoc. 



K 5/3 (k)dk, 



(10) 

(11) 
(12) 
(13) 



where e e is the fraction of energy carried by the shocked 
electrons to proton energy, p is electron energy distribu- 
tion index (see discussion below) which is typically about 

2 — 3 and K§/ 3 (k) is the Bessel function. Finally, is the 
electron energy distribution. Its expression is determined 
by the relations hip between 7 fi .m in 1 7c, max and y c , where 
7 C is defined by (pari et al.lll998h 



7c 



6nm e c 
cjtTB'H 



According to iHuang fe Chens! (|2003h . the expression of 

dN' e 



dN' ■ . , 

-nf- is given by 



Case 1: 7 C < 7e 



die 



= Ci( 7e - l)-^ 1 ), where 



C x = 



P 



(7c,min - 1) P - (7c, me 

Case 2: 7 c , m in < 7c < 7c ,1 



1)~P 



dN' e \ C 2 {le - l)- p for 7o , min < 7e < 
d-fe I C , 3 (7c - l)-(f +1 ) for 7 c < 7c < 7c, n 



(15) 
(16) 

(17) 



Ci 



where 

'(7c,min-l) 1 - p -(7c-l) 1 - p 



+ 



(p-l)( 7c -l) 
(7c - 1) _P - (7c, max 



1)' 



Co 



7c- 1' 
Case 3: 7 C > 7 ,max, 



dN * r 

— = C 4 (7e 

die 



C 4 



where 



p-1 



(7c,min - l) 1 P - (7c, r 



-N K . 



(18) 



(19) 



(20) 



(21) 



In order to calculate the flux at frequency v in the 
observers frame, we first assume that the radiation is 
isotropic in the comoving frame, i.e. 



dP'jv 1 ) 
dVL' 



47T 



(22) 



Let O be the angle between the velocity of the emitting 
m aterial and the line of s ight a nd define /i = cosG. Then, 
bv lRvbicki fc Lightmanl (|l979h . 



1 



dP{i>) _ 
dQ ~ r 3 (l - /3^) 3 4tt ' 

and the relationship between v and v 1 is given by 



(23) 



(24) 



r(i - pn) 

Therefore, the observed flux at frequency v from a source 



with luminosity distance Dl is given by ([Huang et al 
2000a) 



s v = 



i P'(r(i-^) 

AirD 2 L r 3 (l-/3 M ) 3 ' 



(25) 



Although this mechanism is similar to the one respon- 
sible for the afterglow of a GRB, the initial conditions of 
the jet emitted during the accretion to the supermassive 
black hole at the galacti c centre can be very diffe rent from 

" (120061) and iLu et al 



DiacK noic at tne galactic centre can p e very 
(14) that o f a typical GRB. ICheng et al.l {2006) 



(2006) lave estimated the maximum energy AE p carried 
away by the relativistic protons in the jet from the accre- 
tion disk as 



AE P - 6 x lO 52 (77 P /lO- 1 )(M,/M )erg 



(26) 



where rj p is the conversion efficiency from accretion power 
into the energy of jet motion and M* is the mass of the 
disrupted star. ^From this equation, we assume that the 
ejected energy is of the order of about 10 52 ergs. This esti- 
mation of energy in the jet is consistent with the observed 
high energy data in the galactic centre. In case that r) p 
is much smaller than 0.1, we may still be able to observe 
the jet emission if a more massive star is captured. For 
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example, iLu et al. (2006) have studied the situation when 
a red giant star is captured by the black hole. 

The initial Lorentz factor T of the jets ejected from the 
galactic black holes should also be much lower than the 
jets responsible for the GRB afterglows. The jet's initial 
Lorentz factor is > 100 for a typical GRB. However, for a 
typical AGN, the jet's initial Lorent z factor is onl y of th e 
order of several. For example, from ICheng et all (jl999l ). 
the Doppler factors of the studied AGNs are all found 
to be less than 10. For small viewing angle, the Doppler 
factors and the Lorentz factor s of the jets should be of the 
same order. lYuan et al" ( 2002 ) has also estimated the early 
bulk Lorentz factor of NGC 4258 (a low luminosity AGN) 
is about 3. Since during the flare emission, the situation in 
the galactic centre may be similar to an AGN, we should 
restrict the jet initial Lorentz factor to be of the order of 
several. 

4.3. The Results of Numerical Simulation 

We give the numerical simulation of the jet model here 
for illustrating the compatibility between the jet model 
and the flare events. We have used the jet program ini- 
tially for simulating the afterglow light curve of GRB and 
the det ails of the program can be fou nd in iHuang et al.1 
(|2000al| bh and IHuang fc Chengl (|2003h . However, to suit 
the environment in the galactic centre, we have used the 
input parameters as follow: The most important input pa- 
rameters are the jet energy Eoj et = 10 52 ergs, the initial 
bulk Lorentz factor T = 2.5 and the initial jet opening 
angle 9 — 0.3. Other parameters include, the electron en- 
ergy distribution index p — 2.2, the fraction es = 0.01 
of magnetic energy density to thermal energy density, the 
fraction e e = 0.1 of energy carried by the shocked electrons 
to proton energy and the environmental density n = 1000 
cm~ 3 . These parameters are used for all the simulation 
results below unless otherwise sta! ted. Note that they in 
principle can be time and radius dependent, we assume 
that they are constant for simplicity in illustration. We 
also assume that we are looking at the jet along the cen- 
tral axis and the source is 50 Mpc away. In this simplified 
model, we also assumed that all the ejecta are ejected at 
once and there is only one ejection, while in many GRB af- 
terglow model, there can be continuous energy injections 
or several discrete energy injections. We calculated the 
flux in the energy range of 0.1 - 2.4 keV for comparison 
of the results in literatures. 

Our theoretical light curves are shown in Figure [1] 
From our simulation results, when p = 2.2, we found that 
the light curve follows the relationship L x oc t~ a , where 
a changes from ~ 2 at the beginning of the flare to about 
~ 1.3 at the later stage (i.e. about a few years after the 
burst). We can change the power law index a to a broad 
range easily by varying the electron distribution index p. 
Observed data from NGC 5905 are also included for com- 
parison. Note that the intrinsic trigger time of this event 
is unknown. In Figure[T^, we assume that the trigger time 



is about 50 days before th e first observati o n, wh ich is the 
trigger time suggested by lHalpern et al.1 (|2004l ). We see 



that the dash-dotted line can fit the observations well. In 
Figure [TJd, the trigger time is assumed to be about 160 
days before the first observation. Again we see that the 
observations can be well fitted by changing the model pa- 
rameters acco! rdingly. As a conclusion, we found that 
even the jet model alone can explain the observed light 
curve of NGC 5905 reasonably well in both case. 

Furthermore, the above results also illustrate that the 
jet model relaxes the tight constraint on when the flare 
begin since Figure [T^ and [Tp together show that different 
set of model parameters allows the starting time of the 
flare to shift by more than a hundred days. The flexibility 
of the jet model comes from the fact that the realistic light 
curves of the flares actually depend on many parameters, 
such as p and e e , these parameters can pick a range of 
values and it is possible that these values are function of 
time t and position R. Note that if the actual trigger time 
is changed, it mainly changes the time position of the first 
data point in Figure [IJ but it does not affect the positions 
of the other data points significantly. 

In our model, the emission size is much larger than an 
accretion disk. Time evolution of the size of the emission 
regions is shown in Figure [5J However, by the angular 
resolution of the current X-ray telescope, emission size 
is still non-resolvable. Therefore the current observational 
results cannot be used to distinguish if our model is correct 
as long as the source is point-like. 

Figure [3] shows a typical broad band spectrum 
of synchrotron radia tion as described in, for example, 
Granot &: Saril (J2002J) . The synchrotron spectrum is a bro- 
ken power law spectrum with the photon index in the X- 
ray band given by — (p/2 + 1), where p is the energy distri- 
bution index of electron inside the jet. If we take p = 2.2, 
the photon index should be about —2.1, which matches 
spectrum index at the late stage. This gives evidence that 
at least part of the radiation from the flare is from the jet 
rather than solely from the accretion disk. 



5. Discussion and Conclusion 

As a conclusion, in this paper we have summarized the ob- 
served features of some tidal eruption candidates and con- 
cluded that the pure accretion disk model cannot explain 
the radiation properties from these sources completely, es- 
pecially at the late stage of emission. In order to explain 
the situation more throughly, we have proposed the jet 
model. We suggest that when a star is captured by a su- 
permassive black hole at the galactic centre, the star will 
be accreted to the black hole and a jet will be emitted 
during the accretion process. The jet will then decelerate 
in the ISM, and synchrotron radiation is given out in the 
shock emission region. This form part of the radiation in 
the X-ray transient event. This jet model is able to explain 
the late non-thermal spectrum and also relaxes the tight 
constraints on the trigger times of the events. 
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Fig. 1. (a) Observed (0.1 - 2.4 keV) X-ray flux in 
NGC 5905 together with the simulatio n results. The dat a 
points for NGC 5905 is extracted from I Bade et al. ( 19961 ): 
Halpern et al.l ( 2004 ) and the trigger time is assumed to 
be ~ 50 days before the first observation; In our theo- 
retical calculation, we have taken Bo,jet = 3 x 10 51 ergs, 
T=5, e B = 0.001 and e e = 0.06. (b) Same as (a), but the 
trigger time is assumed to be ~ 160 days before the first 
observation, with the simulated light curve is given by an- 
other set of model parameters (i?o,jet = 10 52 ergs, T = 8.3, 
e B = 0.08, e e = 0.5 and p = 2.4). This illustrates that dif- 
ferent parameter sets allow for different trigger time. 



In our calculations, we take the model parameters by 
referring to the situation in AGNs, as discussed in Section 
4.2, since the parameter values are unknown for these flar- 
ing galaxies. In other words, we assume that during the 
flaring states, these "normal galaxies" are similar to the 
AGNs. This is the best choice that we can make currently. 
In the future, when more is known for these flaring galax- 
ies, the model parameters may be refined accordingly. Our 
simulation results show that the jet model alone can ex- 
plain the observed light curve of NGC 5905 reasonably 
well. The jet model also predicts larger emission region, 
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Fig. 2. Time evolution of emission size in different ISM 
density [in cm~ 3 ] 
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Fig. 3. Broad band spectrum predicted by the jet model. 
The used parameters are Eoj e t = 10 52 ergs, T — 5, 9q = 
0.3, p = 2.2, e B = 0.0001, e e = 0.1, n = 1000 cm" 3 . 



though current X-ray telescopes do not have good enough 
angular resolution to test our prediction. 

Although the model predicts emission in all energy 
bands, we may not be able to detect radiation in the op- 
tical band. The main reason for the non-detection may be 
the strong optical background radiation. Note that in the 
jet model, the jet is emitted at the galactic centre, which 
is usually optically bright. Again, NGC 5905 is used as 
the example for illustration (sec Figure QJ. This is the 
only source with optical observation close to the detected 
peak of the flare. Komossa fc Bade] ( 19991 ) concluded that 
during the flare, the optical magnitude of NGC 5905 does 
not varies by more than 0.2 magnitude. This is illustrated 
by the two dashed lines in Figure |4] which has plotted 
the predicted flux in B band using the same parameters 
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used in Figure QJi. For this particular set of model pa- 
ra meters, the model predic tion is consistent with Fig 2 
in K omossa fc Bade! |l999l) as long as the trigger time 



at least 10 days between two consecutive optical obser- 
vations. This example! illustrates that for galaxies with 
optical observations close to the start of the flare, the opti- 
cal observations may help in constraining the trigger time. 
For the other flaring sources, the situation was similar ex- 
cept that the optical observations were done years after 
the peak of the flares, therefore it is even more difficult 
for us to observe the increase in optical flux. 




10 100 
Timefdays] 



1000 



Fig. 4. Optical light curve in B band of NGC5905 pre- 
dicted by the jet model, using the same parameters used 
in Figure la. The Solid horizontal line represents the 
observed optical background (Data taken from Simbad, 
catalog VII/1B) and the dashed horizontal lines repre- 
sent the allowed variati ons based on the results given in 
Komossa fe Bade! (|l999h . 



Another reason may be related to the dust extinction 
There is evidence the dust extinction is v ery serious in 



galac tic nuclei, as in our Milky Way (e.g. iScoville et al 
20031) . In general, dust extinction is important in both 



the optical and x-ray band, however, the way which the 
dust affects the optical and x-ray band is not the same. 
According to IWeingartner fc Draine (|200ll ) , there are evi- 
dence that different galaxies may have different grain sizes 
and it is found that galaxies may have large grain sizes if 
the optical extinction is strong. In the X-ray band, the 
grain size is not so important to the degree of dust ex- 
tinction. Therefore, it may be the case that the flaring 
galaxies are of large grain size if the non-detection of the 
optical band is due to dust extinction. Moreover, if the 
flaring galaxy has a gas/dust ratio that is heavily sub- 
Galactic, then the optical emission will also be seriously 
obscured while X-rays are less affected, since X-ray pho- 
tons are mainly absorbed by gas. 

In the radio band, we find that emission is very weak in 
general (c.f. Fig. [3]). We will not have any detection in the 



radio band unless the source is very near us. However, fu- 
ture radio telescope (e.g. Atacama Large Millimeter Array 
(ALMA)) may be able to detect the size of the emis- 
sion region, by which we can differentiate the region from 
the disk or from the shock emission region. The sources 
are also non-detected in the gamma-ray band, but this is 
simply because the predicted gamma-ray radiation is too 
weak for current detection unless they are observed at the 
very early part of the flare. The thresholds of EGRET and 
GLAS T in cluded in Figure [ 3] were taken from lCheng et all 
(|2004l) and lMcEnerv et alj J2004 . With the high sensitiv- 
ity of GLAST, which will be launched at the end of 2007, 
we should be able to detect the increase in gamma-ray 
flux. 

We noted also that, for most of the identified flares, 
the spectrum during the flare is extremely soft, which 
is not a consequence of synchrotron radiation. In these 
cases, we believe that the synchrotron radiation in the 
jet is still involved in the radiation, but may not be the 
dominated mechanism initially. However, as times goes 
on, the spectrum of the flares become hardened and the 
photon indexes matched with the typical values given in 
synchrotron radiation. By this time, synchrotron radiation 
may become dominant. 

Our jet model allows rebrightening or fast variability 
to occur in the X-ray light curve. For example, when there 
is a density jump in the ISM or when there is a second 
injection of energy or materials into the jet, there could 
be a sudden relati vely small scale ri se in luminosity in 
the light curve (e.g. Huang et al. 20061 ). However, it is un- 
lucky that the current observational data are only sparsely 
sampled so that these fast variabilities have not been re- 
vealed. Moreover, there could be many different ways that 
a jet can evolve in the ISM. As an example, the geomet- 
rical shape of the jet may be cylindrical, but not conical 
as assumed in our current study. Actually, cylindrical jets 
have been observed in a few radio galaxies. The emission 
from a conical jet and a cylindrical jet can be very dif- 
ferent ( Cheng et al. 200 ll ). It worths further investigation 
on which kind of jet works better in explaining the X- 
ray flare events. The results can help us understand the 
environment of the galactic centres. 
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